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We summarize recent developments in our understanding of low-mass quark-antiquark 
excitations in hadronic matter under various different conditions. This includes the ther- 
modynamics of the chiral condensate, pions as Goldstone bosons in normal nuclear matter, 
and excursions into extreme territory of the QCD phase diagram: lepton pair production 
from a fireball expanding through the transition boundary between the quark-gluon and 

t^T" hadron phases of QCD. 

O 

1. Phases of QCD 

Exploring the QCD phase diagram is undoubtedly one of the great challenges in the 
physics of strong interactions. At temperatures exceeding T c — Aq C d ~ 0.2 GeV, one 
expects a plasma of quarks and gluons released from their confining forces. At T < Tc 
and at moderate baryon densities, the relevant QCD degrees of freedom are color-singlet 
hadrons. Chiral symmetry is spontaneously broken. The vacuum is a condensate of 
scalar quark-antiquark pairs. Pions act as Goldstone bosons. Their decay constant, 
/„- = 92.4 MeV, determines the chiral scale ^irf^ ~ 1 GeV which governs the low mass 
hadron spectrum. The lightest vector mesons (p, u) can be seen as the lowest qq "dipole" 
excitations of the QCD vacuum. Current algebra combined with QCD sum rules |ij 
connects their masses directly with the chiral scale, \f2my = 4^/^, to leading order (and 
in the large N c limit). 

Evidently, investigating the changes of the spectral distributions of pseudoscalar and 
vector (as well as axial vector) excitations of the QCD vacuum with changing tempera- 
tures and baryon densities, from moderate to extreme, is a key to understanding QCD 
thermodynamics, its phases and symmetry breaking patterns. 

2. Chiral thermodynamics: selected topics 

A central point in the discussion of the QCD phase diagram is the transition from the 
Nambu-Goldstone realization of chiral symmetry (with non-zero condensate (qq)) to the 
"restored" Wigner-Weyl realization in which the chiral condensate vanishes. Lattice QCD 
indicates that chiral restoration is probably linked to the transition between composite 
hadrons and deconfined quarks and gluons. 
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Spontaneous chiral symmetry breaking (together with explicit breaking by the small u- 
and d-quark masses, m u ^ < 10 MeV) implies the PCAC or Gell-Mann, Oakes, Renner 
(GOR) relation, m^/J = —m q {qq), to leading order in the quark mass m q = {m u + m d )/2. 
The GOR relation continues to hold || in matter at finite temperature T < Tq and 
baryon density p, when reduced to a statement about the time component of the axial 
current. For this component we denote the corresponding in-medium pion decay constant 
by /;(p,T). One finds [| 
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where (qq) Pt T now stands for the T- and p-dependent chiral condensate. The "melting" of 
this condensate by heat or compression therefore translates primarily into an in-medium 
change of the (timelike) decay constant of the pion, given that its in-medium mass m* is 
not much different from m n in vacuum, because of its Goldstone boson nature. 

The chiral scale, 471"/* (p, T), is expected to decrease when thermodynamic conditions 
change toward chiral restoration. This scale defines a gap in the low-energy hadron 
spectrum. A decreasing chiral gap should imply characteristic observable changes in the 
meson mass spectrum. 

Fig. |T| shows typical examples of calculated vector meson spectral distributions in nu- 
clear matter at T = and, in contrast, their evolution with temperature at zero baryon 
density. Such spectra are used in calculations of lepton pair production rates from ul- 
trarelativistic heavy-ion collisions as described in later sections. The main part of this 
presentation focuses on the quest for signals of the QCD equation of state through dilep- 
ton radiation from hot matter. Before exploring such extreme territory, we summarize 
two topics of current interest in chiral thermodynamics at more moderate conditions. 

2.1. Thermodynamics of the chiral condensate 

Suppose we start from a chiral effective Lagrangian, appropriate for the hadronic phase 
of QCD, with Goldstone bosons (pions) coupled to baryons (nucleons). Let Z be the 
partition function derived from this theory, and p the baryon chemical potential. The 
pressure is p(p, T) = (T/V)lnZ, where V is the volume. Given this equation of state, a 
variant of the Hellmann-Feynman theorem (with the quark mass treated formally as an 
adiabatic parameter) in combination with the GOR relation leads to the following density 
and temperature dependent chiral condensate M: 
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with the baryon and scalar densities, p = dp/dp and p$ = —dp/dM, and the nucleon 
sigma term a^ = m q dM/dm q ~ 0.05 GeV, where M is the nucleon mass. 

At T = we have m q [(qq) p — (qq)o\ = UnPs- At low baryon density where ps — p, the 
magnitude of the condensate decreases linearly with density. This behaviour persists up 
to about the density of normal nuclear matter. Detailed calculations using self-consistent 
thermal effective field theory at the two-loop level 0Q] show that effects non-linear in p 
take over at higher density such that the rapid decrease of \{qq) p \ is slowed down. 



The temperature dependence of the condensate at p 
fluctuations. The result of our recent calculations SI 



is determined by thermal pion 
is close to that found in chiral 



perturbation theory || and quite similar to the result of lattice QCD ||, with a critical 
temperature Tq — 180 MeV (for two flavours). The tendency towards "chiral restoration" 
indicated by the dropping condensate should have observable consequences, not only at 
extreme densities and temperatures, but also already in normal nuclear systems. 

2.2. Pionic s-waves in the nuclear medium 

Let us have a brief look at Goldstone bosons in the medium. In the exact chiral 
limit with m q = 0, pions are massless, and this feature persists at all temperatures and 
densities as long as one stays within the Nambu-Goldstone phase of spontaneously broken 
chiral symmetry. Explicit symmetry breaking causes characteristic deviations from this 
pattern. They become particularly interesting in highly asymmetric nuclear matter, with 
large excess of neutrons over protons. 

There has been a recent revival of interest in s-wave pion-nucleus interactions, for 
several reasons. First, the observation of deeply bound pionic atom states in Pb isotopes 
J7| has sharpened the quantitative constraints on the s-wave optical potential. Secondly, 
there is renewed interest in the theoretical foundations of this optical potential from the 
point of view of chiral dynamics |||8||J . 

To leading chiral order, at T = and in the low density limit, the self-energy for a ir + 
or 7T _ at low energy u and momentum q = is simply given by the Weinberg- Tomozawa 
theorem: 

n±(o;, q = 0) = ±^( Pp - P n) + -, (3) 

so that the primary medium effect is a splitting of the 7r + and ir~ masses, Am(7r ± ) = 
±(p p — p n )/4f%, in asymmetric nuclear matter. Systematic calculations of the ir + and tc~ 
mass shifts up to two-loop order using in-medium chiral perturbation theory || predict 
Am(7r~) ~ 14 MeV and Am(7r + ) ~ —1 MeV at nuclear matter density and at a neutron- 
to-proton ratio N/Z = 1.5 characteristic of the Pb region. The repulsive shift found for 
the 7T~ mass is however only half of what is required by the data on deeply bound (Is and 
2p) states in pionic Pb |7J. Several options to repair this discrepancy exist. We believe || 
that an appealing one is the reduction of the pion decay constant in medium, following eq. 
(1). The replacement /„- — > f* in the pion self-energy comes naturally because in-medium 
chiral perturbation theory must now be performed relative to a modified vacuum with 
shifted chiral condensate. Systematic investigations of deeply bound pionic atoms with 
various isotopic chains of heavy nuclei will provide further insights into these challenging 
questions. 

3. Dilepton production rates and thermal QCD 

Lattice simulations with three light quark flavours indicate that QCD undergoes a 



phase transition at a temperature Tc ~ 150 MeV |T(J. It is hoped that it is possible 
to create the quark-gluon plasma (QGP) phase above Tc in ultrarelativistic heavy-ion 
collisions at CERN and RHIC Dileptons (e + e~ and p + p~ pairs) are interesting probes 
in this context since they do not interact strongly but escape unthermalized from the hot 
and dense region formed in such collisions, the fireball. As the fireball expands, it cools 
off and hadronization sets in at Tc- In the hadronic phase, the main dilepton sources 



are pion and kaon annihilation processes which are enhanced due to the formation of the 
light vector mesons p, uj and <fi. 

The differential dilepton emission rate from a hot domain in thermal equilibrium is 
given by 

dN a 2 Imf[(q,T) 
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where II denotes the spin-averaged time-like self energy of a photon in the heat bath, q is 

2 

the photon four- momentum, a = |-, and the lepton masses have been neglected. Later we 
compare our results with the CERES/NA45 data taken in Pb-Au collisions at 160 AGeV 
(corresponding to a cm. energy of ~ 20 AGeV). This is done by integrating eq.(^j) over 
the space-time history of the collision and, taking into account the experimental detector 
acceptance, over the transverse momentum px and average over the rapidity rj. 

3.1. Quasiparticle description of the quark-gluon phase 

Thermal QCD perturbation theory is presumably not applicable to evaluate ImIT, even 
at high temperatures. A way to proceed is to use input from finite temperature lattice 
simulations. We have shown recently [TT| that it is possible to describe the equation of 
state (EoS) of such systems to a very good approximation in terms of a non-interacting gas 
of quasiparticles with thermally generated masses, incorporating confinement phenomeno- 
logically by a temperature-dependent effective number of active degrees of freedom. 

From asymptotic freedom, we expect that at extremely high temperatures the plasma 
consists of quasifree quarks and gluons. Perturbative calculations find spectral functions 
of the form 8{uj 2 — k 2 — m 2 (T)) with m(T) ~ g s T. As long as the spectral functions of 
the thermal excitations at lower temperatures still resemble qualitatively this asymptotic 
form, a quasiparticle description is expected to be applicable. The thermal excitations 
can then be described by a dispersion relation uj 2 (k) = k 2 + m 2 (T). 

We assume that the thermal quark and gluon quasiparticle masses still behave as g(T)T, 
with an effective coupling strength parametrized as 

g(T) ~ g (l - ^V . (5) 



T 

This form is guided by lattice results which indicate that the phase transition is weakly first 
order or second order, suggesting approximately m ~ (T—Tc) 1 with some "pseudocritical" 
exponent 7. 

The second important element in this model is a temperature-dependent confinement 
factor which reduces the number of thermally active degrees of freedom as T approaches 
To- In practice, the quasiparticle parameters are determined by reproducing lattice results 
for the entropy density. The pressure and the energy density then follow accordingly and 
are in excellent agreement with lattice data whenever a comparison can be made JIT] . 



In the quark-gluon phase, the basic mechanism that produces e + e~ pairs is qq annihila- 
tion. Assume now that these quarks and antiquarks are thermal quasiparticles, and that 
these quasiparticles couple to photons in the same way as bare quarks. We can then use 
the standard one-loop result for ImYi, with bare quark masses replaced by the thermal 
masses, m q {T). The essential QCD dynamics is supposed to be incorporated in these 
quasiparticle masses as they are in accordance with the lattice QCD equation of state. 



3.2. The hadronic phase 

Below Tq, confinement sets in and the effective degrees of freedom change to colour 
singlet, bound qq or qqq states. The photon couples now to the lowest-lying 'dipole' 
excitations of the vacuum, the hadronic J p = 1~ states: the p, u and </> mesons and 
multi-pion states carrying the same quantum numbers. The electromagnetic current- 
current correlation function can be connected to the currents generated by these mesons 
using an effective Lagrangian which approximates the SU(3) flavour sector of QCD at 
low energies. We use the improved Vector Meson Dominance model combined with chiral 
dynamics of pions and kaons as described in [12| . Within this model, the following relation 
between the imaginary part of the irreducible photon self-energy Imfi and the vector 
meson self-energies Hy{q) in vacuum is derived: 



v 9v 



(6) 



Here gy is the ^V coupling constant and Fy is the form factor for the coupling of the 
vector meson to the multipion and KK continuum. 

Finite temperature modifications of the vector meson self-energies appearing in eq.(^]) 
are calculated using thermal Feynman rules. The explicit calculations for the p- and 4>- 
meson can be found in ref.[|l3j. At the one- loop level, the p and <fi are only marginally 
affected by temperature even close to Tc because of the comparably large pion and kaon 
masses: m n ~ Tc, vhk — 3 Tc- The thermal spectral function of the w-meson has been 
discussed in detail in |TJ|]. Here, the reaction ww — > nn was found to cause a considerable 
broadening of the u spectral function, leading to T a; (T C ') ~ 7 r w (0). The resulting photon 
spectral function is displayed in figure [I] (left panel). There is still considerable stopping 
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Figure 1. The photon spectral function R(u) = —12tt/uj 2 ImII(a;) at finite temperature 
and p = (left panel), and at T = and finite density (right panel). The qq line in the 
left panel shows the continuum spectral function in the quark-gluon phase for massless 
quarks. 



of the interpenetrating nuclei at SPS energies, with a net baryon density ps in the central 
rapidity region. At RHIC, on the other hand, first measurements indicate that finite 
baryon density effects should not play an important role. For the evaluation of density 
effects at SPS, we use the results discussed in |I5| . The photon spectral function at finite 
density and zero temperature is depicted in figure (right panel). 

3.3. Fireball model 

The fireball model is set up as follows: We assume a spatially averaged time and density 
profile throughout the evolution of the fireball at any given timeslice. Furthermore, we 
require the evolution of the system after the very initial stages to be isentropic. Given 
the value of the total entropy 5*o and the volume expansion of the fireball V(t), the 
temperature profile T(t) can then be constructed via the EoS of the system, depending 
on the relevant degrees of freedom. For temperatures above Tq<, we use the results of the 
quasiparticle model described earlier. In the region below Tc, an interacting hadronic 
gas is an adequate description of the system. We parametrize our insufficient knowledge 
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Figure 2. Left panel: The temperature dependence of the entropy density s in the 
RHIC scenario as compared to the ideal quark-gluon gas and the noninteracting massless 
pion gas limit (dotted). The three relevant regions used in the model calculation are 
given as ideal hadronic resonance gas (full), interpolation (dash-dot) and quasiparticle 
model (dashed). Right panel: Temperature dependence of the entropy density for SPS 
conditions, including the interpolation to the fugacity corrected value at freeze-out. 



close to Tq by interpolating smoothly between two regimes. This approach is supported 
by the general shape of the EoS on the lattice for two light and one heavy quark, where a 
smooth crossover is indicated [|1(J . The entropy density in our model is shown in figure 0. 
Additionally, at SPS conditions, finite chemical potentials for mesons and baryons have 
to be taken into account. The evolution is stopped in the model as soon as a common 
freeze-out temperature Tf ~ 125 MeV is reached. Given this basic framework, we adjust 
the remaining parameters of the model in such a way as to reproduce several hadronic 
observables such as the rapidity distribution of the produced hadrons, p 4 -spectra and 
particle ratios (for details see HT 



3.4. Results 

Initial temperatures are quite high in our approach (250 MeV for SPS conditions, 300 
MeV for RHIC), and we find a significantly prolonged lifetime (10 fm/c at SPS and 15 
fm/c at RHIC) of the QGP evolution phase of the fireball as compared to the results 
obtained by other groups fll8| . This is a consequence of using the more realistic EoS of 
the thermal quasiparticle model rather than the one of the ideal quark-gluon gas. 
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Figure 3. Dilepton invariant mass spectra for the SPS and RHIC conditions realized in the 
CERES/NA45 experiment |19|] (data, upper panel) and the PHENIX detector (prediction, 
lower panel). 



The results for the fully integrated dilepton rates are shown in figure [|. Comparing to 
the CERES experiment, the model achieves a good overall description. We are also able 
to describe the existing data in separate pj-regions (larger or smaller than 500 MeV). 

In figure £|, the different stages of the fireball evolution building up the dilepton yield 
are resolved in small time steps. For early times, the only contribution comes from the 
qq quasiparticle anihilation processes out of the quark-gluon phase. The movement of the 
invariant mass threshold reflects the temperature dependence of the quasiparticle mass, 
which gets light near the phase transition at t ~ 10 fm. One observes that, in spite of the 
growing fireball volume, the contributions from later timeslices to the total QGP yield 
become successively less important. This behaviour is enforced by the confinement factor 
which reduces the thermodynamically active degrees of freedom near the phase boundary. 
For times later than 10 fm/c, the system enters the hadronic evolution phase without 
going through a mixed phase. The most prominent feature is the rapid filling of the low 
invariant mass region through the density-broadened p meson. 



4. Summary 

In our journey through different regions of the QCD phase diagram, we have first 
discussed the temperature and density dependence of the order parameter of spontaneous 
chiral symmetry breaking, the quark condensate. The in-medium evolution of this order 
parameter should manifest itself in a corresponding dropping of the pion decay constant. 
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Figure 4. Time evolution of the (integrated) dilepton yield for the SPS scenario, as a 
function of dilepton invariant mass fjj 



This in turn may lead to observable consequences for low-energy pion-nucleus interactions 
in systems with large neutron excess. New high-precision experiments producing deeply 
bound states of pionic atoms are a promising source of information to investigate these 
issues. 

The major part of this presentation concentrated on the more extreme conditions en- 
countered in heavy-ion collisions at CERN and RHIC Signals of the expansion from a 
transient quark-gluon phase to the hadronic final states can in principle be observed by 
detecting lepton pairs from the expanding fireball produced in such collisions. We have 
emphasized the importance of using input and constraints from lattice QCD thermody- 
namics in the analysis of such signals. 

We thank N. Kaiser, R. Rapp and J. Wambach for stimulating discussions. 
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